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Abstract 


A  high  temperature  expansion  for  the  partition  function  of  the 
Ising  model  ferromagnet  is  developed  and  the  relations  between  it  and  the 
well  known  low  temperature  expansion  are  found.  Making  use  of  these 
relations  and  various  symmetries  and  special  properties  of  the  coefficients, 
it  is  possible  to  calculate  more  coefficients  in  both  series  and  to  check 
those  coefficients  obtained  by  direct  counting.  The  method  is  applied  to 
the  square  lattice  and  to  cubic  lattices  of  two,  three  and  four  layers. 

The  one  layer  Ising  lattice  is  equivalent  to  a  model  of  mono- 
layer  adsorption  with  nearest  neighbor  interactions,  the  spontaneous 
magnetization  of  the  ferromagnet  corresponding  to  critical  adsorption* 
Complete  and  accurate  adsorption  isotherms  are  obtained  from  the  series 
expansions,  and  qualitative  agreement  with  experiment  is  displayed. 

For  the  many  layer  Ising  lattice  the  approximate  variation 
of  critical  temperature  with  number  of  layers  is  found.  It  seems  likely 
that  the  singularity  in  the  susceptibility  of  the  two  layer  lattice  is 
of  the  same  form  as  that  in  the  plane  square  lattice,  rather  than  being 
intermediate  between  that  of  the  plane  square  and  simple  cubic. 

Certain  mathematical  problems  arising  out  of  this  work  are 


pointed  out. 
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I  INTRODUCTION 


1 .  Cooperative  Phenomena  and  the  Ising  Model. 

The  physical  systems  which  are  treated  by  statistical  mechanics 
may  be  divided  into  two  classes;  those  in  which  the  macroscopic  properties 
are  determined  primarily  by  the  properties  of  the  individual  microscopic 
subsystems  and  in  which  the  constituent  subsystems  may  be  treated  as  being 
essentially  independent,  and  those  in  which  the  interactions  between  the 
subsystems  are  predominant  in  determining  the  macroscopic  behavior.  The 
latter  class  includes  all  systems  which  exhibit  phase  transitions  (with 
the  exception  of  the  condensation  of  a  free  Bose  gas)  such  as  the  conden¬ 
sation  of  real  gases,  ferromagnetism,  order-disorder  transitions,  and 
critical  adsorption.  Since  the  interactions  give  rise  to  correlations 
among  and  cooperative  motion  of  large  numbers  of  the  microscopic  sub¬ 
systems,  the  above  mentioned  class  is  said  to  exhibit  "cooperative 
phenomena" . 


Presumably  the  simplest  possible  model  of  a  cooperative  assembly 
would  be  one  in  which; 

(a)  the  microscopic  systems  are  identical  and  possess  only  two 
possible  states, 

(b)  interactions  are  between  nearest  neighbors  only. 

Such  a  model  was  proposed  for  ferromagnetism  by  Ising  [1]  in  1925.  He 
associated  with  each  lattice  point  in  a  crystal  a  spin  coordinate  <y 
which  took  on  two  values,  +  1  for  spin  "up"  and  -  1  for  spin  "down". 
The  interaction  between  nearest  neighbors  (n  -  n)  was  postulated  to  be 
-  J  <TjO\  ,  that  is  -  J  for  parallel  spins  and  +  J  for  antiparellel 
spins.  In  addition,  there  was  an  interaction  with  the  external  magnetic 
field  H  of  -  H  ,  that  is  -  H  for  spin  with  the  field  and  +  H 


for  spin  against  the  field  (we  take  the  magnetic  moment  per  spin  as 
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un 


ity)*  Thus  the  energy  of  an  Ising  ferramagnet  is 


-  J 


I.  'i  "3 


(n-n)  i 


and  the  partition  function  is 


Z  = 


I  «p  kj  X  ^  +  lf  I'i 

(c, )  L  (n-n)  i 


(1) 


(2) 


where  l.  denotes  the  sum  over  all  states  ^  for 

(<Ti^  V*1  °2=:±1  V*1 


lattice  of  B  sites, 


By  suitable  identification  of  variables  the  Ising  model  can  be 
made  a  model  of  many  different  cooperative  systems.  The  transformations 
necessary  to  treat  the  problems  of  the  lattice  gas  and  the  binary  alloy 
are  dealt  with  in  detail  by  Hill  [2]®  Newel land  Montroll  [3],  and 
Temper ly  [1|.]0  It  was  this  correspondence  between  the  Ising  Model  and 
the  lattice  gass  first  established  by  Lee  and  Yang  [5 ]»  which  enabled 
them  to  show  that  the  condensation  of  a  gas  into  a  liquid  phase  could 
be  rigorously  deduced  from  the  partition  function  for  a  particular  single 
model . 

The  necessary  transformation  to  describe  monolayer  adsorption 
will  be  given  in  detail  in  Ch„  III. 

2.  Motivation  for  Studying  Thin  Films. 

The  experimental  evidence  for  cooperative  phenomena  in  thin  films 

is  quite  varied. 

Firsts  the  experiments  of  Jura9  Loeser,  Basford,  and  Harkins  [6]» 


Fisher  and  McMillan  [7j9  and  Machin  and  Ross  [8],  exhibit  critical 
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adsorption  i„e.  an  abrupt  increase  in  coverage  at  some  pressure  as  long  as 
the  temperature  is  below  a  certain  critical  temper ature? and  smooth  isotherms 
for  temperatures  higher  than  critical.  The  similarity  of  these  results  to 
the  approximate  calculations  of  Fowler  [9]  prompted  the  more  accurate 
treatment  of  this  problem  contained  in  Ch.  III. 

Secondly,  various  papers  have  reported  a  shift  of  a  "A  ~  point”* 
with  the  number  of  layers  of  an  adsorbed  gas.  Experimenting  with  He 

adsorbed  on  FeJ3  ,  Frederikse  [11]  found  a  specific  heat  maximum  at  a 

£  3 

lower  temperature  than  the  bulk  liquid  A  -  temperature,  which  decreased 
in  height  and  position  with  decreasing  coverage. 

Morrison,  Drain,  and  Dugdale  [12]  have  found  a  similar  result 
for  on  Ti  the  transition  occurring  in  a  temperature  range  just 

below  the  bulk  melting  point.  Steele  [13]  has  found  a  somewhat  similar 
effect  in  adsorbed  methane.  These  results  partially  served  to  motivate 
the  study  of  the  many  layer  Ising  model  in  Ch.  If,  which  is  found  to 
exhibit  a  qualitatively  similar  shift  of  A  -  point  with  the  number  of  layers 
although  a  precise  correspondence  between  this  model  and  the  above  mentioned 
systems  does  not  seem  appropriate. 

Thirdly,  the  superfluid  transition  (A  -  point)  in  bulk  liquid 
helium  is  believed  to  be  intimately  related  to  Bose  -  Einstein  condensation. 


*  The  term  "A  =  point"  was  originally  given  to  a  singularity  in  the 

specific  heat  of  liquid  helium,  which  in  early  experiments  appeared 
to  be  a  discontinuous  jump  having  roughly  the  shape  of  a  reversed 
Greek  lambda.  More  recent  experiments  [10]  indicate  that  this  sing- 
ularity  is  really  an  infinity,  and  so  the  term  "A  -  point"  will  be 
used  to  denote  a  singularity  of  this  type. 
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Long  and  Meyer  [14]  have  found  superflow  in  He  films  as  thin  as  -w  2 
atomic  layers „  the  onset  temperature  for  superflow  decreasing  from  the 
bulk  A  -  point  with  decreasing  film  thickness.  Manchester  and  Buckley 
[15]  have  found  a  peculiar  phenomena  of  temperature  oscillations  in  films 
as  small  as  1  or  2  atomic  layers,  which  they  interpret  as  being  due  to 
the  onset  of  super flow  in  these  thin  films.  However  it  seems  to  be 
theoretically  established  that  a  two  dimensional  ideal  B.-E.  gas  does 
not  condense  as  does  a  three  dimensional,  gas  [16,  17]. 

This  has  given  rise  to  the  conjecture  (not  yet  settled)  that 
the  critical  properties  of  thin  He  films  may  be  due  to  an  order~ 
disorder  transition. 

For  these  diverse  reasons  it  was  considered  interesting  to 
investigate  the  application  of  the  Ising  model  to  thin  films. 

3.  Scope  of  this  Thesis. 

The  material  of  this  thesis  is  in  three  main  sections.  In 
Ch„  II  some  new  work  on  the  general  theory  of  the  Ising  model  is  presented, 
after  reviewing  the  results  of  previous  authors.  Then  in  Ch.  Ill  and  I? 
this  theory  is  applied  to  particular  lattices  of  interest.  The  final 
chapters  contain  a  discussion  of  questions  which  are  raised  by  the  work 
of  this  thesis  and  a  summary  of  conclusions. 
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II  THEORY  OF  THE  IS  INC  MODEL 

4.  Review  of  the  Literature. 

The  very  extensive  literature  on  the  Ising  problem  has  been 
covered  in  two  review  papers  by  Neweliand  Montroll  (1953)  [3]  and  Domb 
(i960)  [18].  Hill's  book  (1956)  [2]  also  contains  a  good  review.  We 
shall,  therefore,  merely  summarize  the  major  results. 

The  model  was  proposed  by  Ising  [1]  in  1925s  who  solved  it  for 
the  one  dimensional  case.  However  this  case  does  not  exhibit  any  ferro- 
magnetic  properties  (spontaneous  magnetization,  curie  point)  at  non-zero 
temperature,  so  attention  was  centered  on  the  much  more  difficult  two 
and  three  dimensional  cases.  Until  1941  the  investigations  were  carried 
out  by  means  of  approximate  methods.  An  excellent  review  of  these  as 
well  as  more  recent  approximate  methods  is  contained  in  ref.  [18].  For 
reasons  given  below  we  believe  that  exact  series  expansions  provide  a 
more  fruitful  approach  than  any  of  these  approximate  methods,  and  so  no 
detailed  references  to  them  will  be  given. 

In  1941  Kramers  and  Wannier  [19]  showed  that  the  partition 
function  could  be  derived  as  the  largest  eigenvalue  of  a  matrix,  and 
for  the  two  dimensional  square  lattice  they  found  a  symmetry  which  per¬ 
mitted  them  to  locate  the  curie  temperature, assuming  its  existence  and 
uniqueness.  This  same  symmetry  also  ruled  out  the  possibility  of  certain 
types  of  singularities  including  a  finite  jump  in  specific  heat  at  the 
curie  temperature,  which  all  of  the  approximate  treatments  had  predicted. 

The  two  dimensional  square  lattice  with  zero  magnetic  field  was 
solved  exactly  by  Onsager  (1944)  [20]  who  found  that  the  singularity  in 
the  specific  heat  is  a  logarithmic  infinity.  A  simpler  but  still  complex 


’ 
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solution  was  given  by  Kaufman  (1949)  [21 ]  making  use  of  spinor  represent¬ 
ations.  As  an  alternative  to  the  matrix  methods  of  Onsager  and  Kaufman 
(OK),  Kac  and  Ward  [22]  (KW)  in  1952  produced  a  solution  from  a  combin¬ 
atorial  method  involving  the  counting  of  closed  graphs.  However  they 
could  not  justify  the  details  of  their  calculation,  and  the  complete 
proof  was  not  achieved  until  i960  by  Sherman  [23].  Recently  (I96I),  the 
Soviet  physicists  Dykhne  and  Rumer  [24]  have  published  an  apparently  much 
simpler  solution  of  the  combinatorial  method  by  taking  a  different  approach 
from  that  of  KW  and  counting  the  graphs  by  means  of  a  formalism  involving 
creation  and  annihilation  operators.*  However  none  of  these  methods  have 
yet  permitted  a  generalization  to  either  non-zero  magnetic  field  or  to 
three  dimensions. 

The  methods  of  OK  were  applied  by  several  workers  independently 
[253  (1950)  to  the  triangular  and  hexagonal  lattices,  which  were  found 
also  to  possess  logarithmic  singularities  in  their  specific  heats.  Yang 
[26]  (1952)  extended  OK"s  method  to  obtain  the  spontaneous  magnetization 
of  the  square  lattice.  This  was  generalized  to  the  triangular  lattice  by 
Potts  [27]  (1952)  and  to  the  hexagonal  lattice  by  Naya  [28]  (1954). 

In  addition  to  the  exact  closed  solutions  which  have  so  far  been 
obtainable  only  for  one  and  two  dimensional  lattices,  it  is  possible  to 
derive  exact  series  expansions  by  methods  which  are  equally  applicable  to 
any  lattice.  Van  der  Waerden  [29]  (1941)  developed  a  high  temperature 
series  for  zero  magnetic  field  involving  the  counting  of  closed  graphs, 
which  was  the  basis  for  the  above  mentioned  combinatorial  method  of 
solution.  His  series  is  generalized  to  arbitrary  magnetic  field  in  this 
thesis.  Series  valid  for  high  and  low  temperatures  and  non=zero  magnetic 
field  were  calculated  by  Domb  [30]  (1949)  for  the  plane  square  lattice 


*  This  was  first  done  by  C. A. Hurst  and  H.S. Green,  J.Chem.Phys. 33.1059(1960) . 
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using  a  complicated  matrix  method.  Similar  series  were  found  for  the 
simple  cubic  lattice  by  Wakefield  [ 31 ]  (l95l)»  who  used  the  simpler  method 
of  configuration  counting.  Series  have  also  been  calculated  for  other 
lattices,  a  listing  being  given  in  the  appendix  of  ref. [32]. 

It  is  possible  to  test  approximate  methods  by  comparison  with 
these  exact  expansions.  This  has  been  done  for  the  plane  square  lattice 
by  Kramers  and  Wannier  [19]  and  for  the  simple  cubic  by  Ter  Haar  [33]. 

In  all  cases  the  approximate  methods  were  valid  to  a  lower  order  of 
accuracy  than  could  be  obtained  from  the  exact  series  without  undue  effort. 
For  this  reason  and  the  fact  that  approximate  methods  are  usually  seriously 
in  error  about  the  nature  and  position  of  the  singularity  we  decided  to 
use  only  exact  methods  in  this  thesis. 

Recent  work  has  been  directed  toward  extending  the  calculation 
of  exact  series.  Domb  and  Sykes  [18,  3^3  have  developed  methods  which 
simplify  the  counting  of  the  graphs  arising  in  these  expansions. 

Rushbrooke  and  Eve  [353  have  made  use  of  a  digital  computor  to  count 
simple  closed  polygons  on  a  cubic  lattice.  When  a  sufficient  number  of 
terms  in  a  series  are  known  it  is  then  possible  to  make  a  reasonable 
estimate  of  the  nature  of  the  singularity  by  extrapolating  ratios  of  the 
coefficients  [36]. 

5.  Series  Expansions.* 

We  shall  here  derive  the  two  series  expansions  of  the  Ising 
partition  function  (p.f)  which  form  the  basis  of  the  calculations  in 
this  thesis. 


* 


The  material  in  this  and  the  next  section  have  been  submitted  for 
publication  in  J.  Math.  Phys.. 


aw  .  ■  . ft* 


-  8  - 


Let 


and 


B  =  no.  of  spins  in  the  lattice  , 

N  =  no.  of  negative  ({)  spins  , 

A  =  no.  of  nearest  neighbor  (n-n)  bonds  in  the  lattice  , 
Y  =  no.  of  antiparallel  (It)  n-n  bonds  . 


Then 


I  ».  ■ 


B  -  2N 


and  ^  cr  ,  u ,  =  A  -  2Y 


(n-n) 


i  j 


Eq.  (2)  becomes,  upon  substitutions 


Z  =  exp  (AK  +  Bl)  ^  exp  (-2KY  -  2IN) 


J  H 

where  K  =  kT  and  1  =  kT 


This  can  also  be  expressed  as 


Z  =  x 


where  x  =  e 


■A/S  - B/2 


•2K 


y  S  e 


B  cN 

I  I 

N-o  Yf=o 

=21 


v  \  Y  N 
,YS»B)  x  y 


(3) 


g(N,YgB)  =  no.  of  arrangements  of  N(|)  spins  forming 
Y(if)  n-n  pairs , 

and  c  =  max.  no.  of  nearest  neighbors  per  site. 

We  shall  denote  the  double  sum  in  (3)  by  the  symbol  (x»y)  • 

-*■  0  as  T  0  we  shall  refer  to  (3)  as  the 


Since  x  and  y 
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low  temperature  (L.T..  )  expansion. 


Equation  (2)  may  also  be  written  as 

Iff. 


n 

n-n 


n  * 

i 


(2°) 


Since  (<r^)2  =  (o^  <t^)2  =  1  s  we  have 

K<TiCT4 

e  =  (1  +  o-^  CTj  u)  cosh  K  s  u  =  tanh  K 

and 

Iai 

e  =  (l  +  v)  cosh  I  ,  v  =  tanh  I  „ 

Thus 

Z  =  (cosh  K)A  (cosh  l)B  ^  (l+o^o^u)  |J  (l+o^v)  .  (i|) 

(cri)n<=n  1 

If  we  identify  a  bond  with  the  function  o\  »  the  product  of  the 

spins  at  its  endss  then  expanding  the  two  products  in  (h)  will  lead  to 

£  m 

a  double  power  series  in  u  and  v  with  the  coefficient  of  u  v 
equal  to  El  (all  combinations  of  i  distict  bonds)  x(all  combine 
ations  of  m  distinct  spins).  A  term  in  this  sum  can  be  represented 
graphically  on  the  lattice  by  drawing  a  line  for  a  bond  and  circling  a 


IS 


spin.  For  example 


1 


. 


10 


represents  the  term  [(ff5°'7)(^7^8)(^7^11)(^10cru)(tr6CTl0)  ]  x  [a^Og]  , 

5  2 

which  occurs  in  the  coefficient  of  u  v 


Now  the  configurational  sum. 


will  yield  a  non-zero 


result  only  for  terms  in  which  every  a  appears  to  an  even  (or  zero) 
power.  Thus  the  contributing  terms  must  be  such  that  the  l  bonds  form 
a  graph  with  m  odd  vertices  and  the  m  spins  coincide  with  these  odd 
verties.  (An  odd  vertex  is  the  meeting  point  of  an  odd  number  of  bonds). 
So  we  obtain 

K  21 

Z  =  (cosh  K)A(cosh  l)B  2B  ^  i**(/,m;B)  vm  u^  (5) 

t—o  m=o 

where  ^(ism;B)  =  no„  of  graphs  of  l  bonds  with  m  odd  vertices 

which  may  be  drawn  on  the  lattice.  We  shall  denote  the  double  sum  in 

B 

(5)  by  the  symbol  Q  (usv)„  The  upper  limit  of  21  on  £  arises 

m 

because  a  graph  of  l  bonds  has  at  most  21  vertices. 


Since  u  and  v  — ►  0  as  T  ->oo  s  (5)  will  be  called  the  high 
temperature  (H.T.)  expansion. 

Equation  (5)  is  a  generalization  to  arbitrary  magnetic  field 
of  Van  der  WaerdenJs  [29]  expansion  for  zero  field,  which  involves  only 
the  closed  graphs.  Oguchi’s  [37]  expression  for  the  zero-field  suscepti¬ 
bility  in  terms  of  graphs  with  two  odd  vertices  is  also  a  special  case 
of  eq.  (5). 


We  can  also  prove,  from  eq.  (5)» 
extensively  by  Domb  [30].  By  substituting 


0 


£\ 

■  1 


p/y) 


I--O 


(i+y) 


2  i 


u 


a  theorem  conjectured  and  used 
v  =  (l-y)/(l+y)  »  we  find 


■ 


•  > 
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where  p^(y)  is  a  polynomial  of  degree  <  21  .  A  proof  was  given  by 
Wakefield  [31]  in  1951  but  eq.  (5)  provides  a  much  simpler  proof. 

6.  Symmetries  and  Special  Properties. 

The  combinatorial  coefficients  in  eq.  (3)  and  (5)  have  many 
useful  properties.  Moreover  since  (3)  and  (5)  represent  the  same 
function,  every  property  of  one  set  of  coefficients  must  imply  a  pro¬ 
perty  of  the  other.  Some  of  these  corresponding  properties  are  collected 
in  table  I.  From  any  one  of  the  properties  in  a  row,  one  can  deduce  the 
other  two  by  equations  (3)  and  (5). 

We  shall  identify  the  columns  of  this  table  by  (p.f.),  (L.T.), 
and  (H.T.),  and  the  rows  by  letters  (A)„.„(e)„ 

Row  (A)  states  that  the  invariance  of  the  p.f.  under  magnetic 
field  reversal  follows  both  from  the  fact  that  reversing  all  spins  does 
not  change  the  number  of  antiparallel  pairs  (L.T.,A),  and  from  the  fact 
that  graphs  with  an  odd  number  of  odd  vertices  do  not  exist  (H.T„,A). 

Row  (B)  states  that  the  thermodynamic  equivalence  of  ferro¬ 
magnetic  and  antiferromagnetic  loose-paeked  lattices  follows  both  from 
the  fact  that  reversing  all  spins  in  one  sub= lattice  changes  (ft)  pairs 
into  (If)  pairs  and  vice  versa  (L.T.,B),  and  the  fact  that  closed  graphs 
of  odd  order  cannot  be  drawn  on  such  a  lattice  (H.T.,B). 

Rows  (C),  (D)  and  (E)  express  the  properties  which  make  simple 
the  limiting  cases  of  zero  interaction,  infinite  magnetic  field,  and 
infinite  interaction  respectively.  The  remaining  limiting  case  of  zero 
magnetic  field,  of  course,  does  not  yield  a  simple  solution. 


The  property  (H.T„,E)  deserves  special  mention  since  it  is  the 


.  .•  -  -  ■ 


■ 


- 


11.1 


TABLE  I  -  CORRESPONDING  PROPERTIES 


Property  of: 


Partition  function 

L.T.  coefficients 

H.T.  coefficients 

Z(K,I)  =  Z(K,-I) 

g(N,Y;B)  =  g(B-N,Y,B) 

Jj/(i,m;B)  =  0  for  m  odd 

(A) 

*  Z(K,0)  =  Z(-K,0) 

B 

V  g(N,Y;B)  3  G(Y;B) 

N^o  •  G(A-YiB) 

J^(J, 0;B)  =  0  for  i  odd 

(B) 

aB  (i»y)  =  (i+y)B 

cN 

Y  g(N.Y;B)  =  (B) 

Y=o 

J 

(0,m;B)  =  1,  m=0 
=  0,  m>0 

(c) 

A8  (*,0)  =  1 

g(0,Y;B)  =  1,  y=0 
=  0,  y>0 

2i 

y  m;B)  =  (A) 

m=0 

(D) 

aB  (0,y)  =  l+yB 

i 

g(N,0;B)  =  1,  N=0  or  B 
=  0  otherwise 

A 

y Mi, mjB)  =  2A'B  (®) 

i=o  (u(-i)m) 

(E) 

*  This  property  holds  only  for  a  loose-packed  lattice,  which  can  be 

divided  into  two  sub- lattices  such  that  members  of  a  nearest 
neighbor  pair  are  always  in  opposite  sub-lattices. 
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only  property  in  this  table  which  could  not  be  seen  by  direct  inspection 
and  was  in  fact  discovered  by  inference  from  the  corresponding  L.T.  property. 

It  is  convenient,  and  eventually  necessary,  to  pass  to  the  limit 
of  an  infinite  lattice.  From  table  I  (L.T.,C)  and  the  fact  that  all 
g(N,Y;B)  are  positive  Increasing  functions  of  B,  we  see  that  for 
0  <  x,  y  <  1  , 

g 

A  (x,y)  is  an  increasing  function  of  B  and 

AB(x,y)  <  AB(1,1)  =  2B  . 

From  eq.  (3) 

£inZ«-||inx--iny  +  g-in  AB(x,y) 
£-||inx-|iny+in2  . 

A 

Thus  for  any  reasonable  lattice  in  which  7  *  jjlim  3  exists  we  will  have 

lim  ^  in  Z  exists.  This  can  also  be  proven  from  eq.  (5)  and  table  I 
B-+00  B 

(H.T. ,D) ,  preserving  the  symmetry  between  (3)  and  (5). 


Taking  this  limit  in  (3)  and  (5)  leads  to 

lim  -  in  Z  b  -  -  7  in  x  -  -  in  y  +  in  A  (x,y)  (3’) 

B  00 

s  7  in  (cosh  K)  +  in(cosh  X)  +  in^  +  infl(u,v) 

(5') 

Upon  subtraction  and  a  little  algebra  we  obtain 

in  A  -  in(l+y)  =  in  fi  *  7  in(l+u)  }  (6) 


an  equation  to  which  we  will  return  later. 


Since  g(0,0;B)  ■  ^(0,0;B)  =  1  ,  the  following  formal 


I  I 


■ 


.  .  3*»«n  v  t  f 


■ 


. 


r 
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expansions  can  be  made: 


In  A  s 


B  CN 

5  <n  {X  Z  *V  ]•. 


N=o  Y=o 


00  cN 


=  Z  Z  ®<n,y)  xYyN  ■ 


N=1  Y=o 


(7a) 


A  2 1 


In  Q  =  -  In  ^  JL,  (/,m;B)  vV*  j- 


l=o  m=o 
00  21 


=  ^  ^  vV 


l-l  m=o 


(7b) 


Now  in  the  equation 


In 


(i+  z 

n=l 


a 

n 


•“>  -l 

n=  1 


b 

n 


n 

z 


it  is  shown  in  appendix  I  that  a ^  and  b^  are  related  by  equations 
of  the  form 


b 

n 


n 


°Sl  -  3-  “cl  o 

n-1  1  n-2  2 


2  2 

+  ar=2al  +ar-4a2  + 


The  significant  properties  of  this  relation  are: 

(a)  If  the  subscripts  are  considered  to  measure  "degree",  then 
all  terms  are  of  "degree"  n  , 

(b)  The  term  aR  is  the  only  one  linear  in  a's  and  its 
coefficient  is  always  unity. 

In  (7a)  it  is  appropriate  to  identify 


i  ■  ... 
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cN 

g(N,Y;B)  xY  with  a 

n 

Y=o 

cN 

and  B  g(N,Y)  xY  with  bn  , 

Y=o 

and  similarily  for  (7b).  Property  (a)  above  justifies  keeping  the  same 
limits  on  the  two  sums. 


Now  if  the  lattice  boundaries  are  cyclically  connected  and 

B  »  N,  A  »  i  (so  that  configurations  which  could  loop  the  lattice 

and  intersect  themselves  do  not  occur)  then  g(N,Y;B)  and  ^(/,m;B)  are 

polynomials  in  B  with  no  constant  term  (N,//o).  Furthermore  we  know 

lim  ^  In  Z  exists.  But  this  will  be  possible  only  if  all  powers  of 
B  — ►  00  B 

B  greater  than  the  first  cancel  out  in  the  transformation  from  the  a’s 
to  the  b's.  This  cancellation  provides  a  check  on  the  consistency  of 
the  coefficients  g(N,Y;B)  and  (i,m;B)  but  it  is  usually  a  trivial 
one. 


Since  g(N,Y;B)  and  m;B)  contain  no  constant  term 

(N,//o),  property  (b)  above  implies 


g(N,Y)  0  coef.  of  B  in  S(N,Y;B) 
and  -  coef.  of  B  In  m;B) 


(8) 


provided  we  have  cyclic  boundary  conditions  and  B  »  N,  A  »  l  . 

I 

Domb  [58]  also  gives  a  proof  of  (8)  but  he  fails  to  state  the  conditions 
for  its  validity  which,  we  shall  see,  can  be  important. 

From  eq.  (8)  certain  useful  results  follow.  The  most  obvious 


is; 


■  -  f' 


. 

3<a  00  u““ 

*>  3  L'?:  'floo  <Ij  ...  ;  '  x:>  f.  '>iXo.D  Of!  >■..•■•  I  c;:  ■  3ti  *  r*K 


y.  '  W-  {  <  :•  :  0  ■  ■■  .  i  :  ■  ;k  • uit;; 


' 


' 
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■  0  then  (i,m)  =>  0  , 

and,  if  g(N,Y;B)  ■  0  then  g(N,Y)  *  0  . 

From  table  I  (H.T.,D)  we  have 

£  »;B)  =  (A)  =  a(a-i)...(a-/+i) 

m=o 


Since  A  ^  7  B  ,  eq.  (8)  yields  (by  taking  coef.  of  B  on  both  sides) 

£  &{t, ■»)  =  1 


J-l 


m=o 

Similarily 

CN 


Y  s(H.Y)  . 

N 

Y=o 


We  run  into  difficulty,  however,  when  we  attempt  to  perform 
such  a  calculation  on  table  I  (H.T.,E)„  For  example  we  consider  m  =  0. 

V  i,U,0;B)  =  2  2A‘B  =  2  e<A-B>/n2 

=  2  (l  +  (A-B)in  2  +  g]  +  ...  1 

Formally  eq.  (8)  yields 
00 

y  jba,  0)  =  2  in  2  (7“l)  in  the  limit  as  A  ->-oo  . 

l—O 

OO 

But  subsequent  calculations  will  make  it  clear  that  diverges. 

£= o 

The  resolution  is  the  fact  that  the  upper  limit  of  the  initial  sum  is  A, 
so  that  the  restriction  A  »  l  can  not  hold  and  (8)  does  not  apply. 


Let  us  now  substitute 


1;  .. 


IS 


•w  'Hktrsvctf  "uii-Vb  c*«J-  *<<  aV 

~j 


-  16  - 


*n  =(ftO“-  I  c”  “r  ■ 


r=o 


Lr.l 

i  +  y 


n 


uu 

■  I 


_  n  r 
C  y 
r  J 


r=o 


(9) 


Into  eq,  (6).  (Equation  (9)  defines  the  expansion  coefficients  „ 

Some  of  their  properties  and  the  method  of  calculating  them  are  contained 
in  appendix  II). 


Collecting  powers  of  u  and  y,  we  obtain 
21  cN 

V  CN"  .  V  g(N,Y)  cj  ,  i  /  0,  N  /  0  , 

m=2  Y-l 

y  ,  l  /  0  , 

m=o 

f  i\N+1 

g(H,Y)  -  -  ,  N  /  0  . 

Y=o 


(10) 


(11) 


(12) 


The  lower  limits  in  (10)  have  been  raised  from  zero  to  Y  =  1  and 
m  =  2  because  Cr°  =  0  for  r  /  0  and  Jb(£9m)  =  0  for  m  odd. 
Equations  (ll)  and  (12)  have  already  been  derived  from  table  I  and 
eq.  (8). 


In  the  following  chapters  the  utility  of  (10)S  (ll),  and 
(12)  will  be  shown  by  application  to  certain  lattices  of  interest. 

7.  Calculation  of  Physical  Quantities  from  the  Partition  Functions. 


It  is  shown  in  the  text  books  on  statistical  mechanics,  e.g. 
Rushbrooke[39],  that  the  equilibrium  value  of  a  quantity  Q  is  given 


■ 


by  the  Boltzmann  average  over  all  microstates  of  the  assembly: 


Q 


?Q1 


e 


-Ej/kT 


E  rEl>T 

1 


Thus  for  the  Ising  Model  we  have,  from  eq.  (l)  and  (2),  the  Internal 
energy  per  spin  given  by 


E 


kT2 


dT 


in  Z 

B 


(13) 


and  the  magnetization  per  spin  given  by 

(Z  a  )  e“Ei AT 


M  = 


Z 

1 


'j 


B  Z 


=  kT 


in  Z 

B 


(lb) 


where  we  have  temporarily  used  the  subscript  i  to  lable  microstates 
rather  than  spins,, 

It  is  necessary  for  calculation  to  express  (13)  and  (14)  in 
terms  of  the  variables  x,  y9  u,  and  v  . 

From  the  definitions  of  these  variables  1 

=2JAT 

x  =  e  '  ,  y  =  e 

u  =  tanh  (J/lcT)  =  ,  v  =  tanh  (H/kT)  =  , 


we  find  [ 3I 3 


and 


6  2J  6 

5t  =  £?2  x  + 


626 

“  ”  kT  y  6y 


2H 

kT2  ^  6y 


(ref.  31  has  omitted  the  minus  sign). 


(15) 

(16) 


Using  (3')  we  then  find  the  zero  field  specific  heat  to  be 


V  :«J£. 
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and  the  magnetization  to  be 


M  =  1  -  2  y  in  A 


(18) 


In  terms  of  H.T„  variables  we  have 


a  Y  i 

x  s;  ;  =  k 


(i 


2\  ^  ”U 

-  u  >  SC  J 


and 


^  (  -t  2\ 

3H  =  kT  ^  ’  V  ^  dv 


Using  (5*)  or  (17)  and  { 18)  with  (b)  we  obtain 


(19) 


(20) 


Cv  =  k  (in  x)2irj^7^y2  +  l  [\l"U2)  J2  in  (21) 


and 


M  -  v  +  (1  =■  v2)  ^  in  fl 


(22) 
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III  MONOLAYER  ADSORPTION  ON  A  SQUARE  LATTICE  SUBSTRATE* 


8.  The  Ising  Model  for  Monolayer  Adsorption. 

The  model  of  adsorption,  which  we  shall  show  to  be  equivalent 

to  the  Ising  model,  is  one  in  which  the  solid  substrate  is  considered  to 

of 

be  a  regular  two  dimensional  array  of  sites,  each  capable^ holding  one 

adsorbed  gas  molecule.  There  is  an  interaction  energy  -E  between  each 

8 

adsorbed  molecule  and  the  surface  and  an  interaction  energy  -V  between 
molecules  adsorbed  on  nearest  neighbor  sites.  Approximate  treatments  of 
this  model  were  first  given  by  Fowler  [9]  and  Peierls  [40],  who  used 
Bethe's  method.  More  recently  Bumble  and  Honig  [41]  have  used  the  rather 
involved  approximate  method  of  Hijmans  and  de  Boer  [b2] .  In  this  treat¬ 
ment  we  shall  show  how  certain  important  features  follow  from  exact 
knowledge  about  the  Ising  partition  function  and  the  complete  isotherms 
can  be  accurately  calculated  from  series  expansions. 

In  order  to  form  the  grand  partition  function  of  the  adsorption 
surface  let  us  assign  to  each  site  a  coordinate  o\  ,  i  =  1  . . .  B  , 
which  takes  the  value  +1  if  the  site  is  occupied  and  -  1  if  the  site 
is  unoccupied.  The  number  of  adsorbed  molecules  will  be 

Na  =  Y  (l  +  <r.)  /  2  ,  (23) 

T 

and  the  number  of  nearest  neighbor  pairs  of  adsorbed  molecules  will  be 

x  =  y  §  +  <r.)  (1  +  <r.)  /  k  ,  {2k) 

n-n 


*  The  material  of  this  chaper  has  been  submitted  for  publication  in 
Physica  in  a  joint  paper  by  L.  E.  Ballentine  and  D.  D.  Betts. 


' 
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- 


. 


* 
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and  the  grand  partition  function  for  the  surface  will  be 

1  ■ 


2  -  E  exp  {  “kT5  + 

(Oj) 


E  N  +  VX 
s  a 


kT 


(25) 


In  expanding  the  expression  for  X  it  is  convenient  to  write 


I  •.  ■  1  E 


3ij  Ti 


n-n 


where  =  1  if  i  and  j  are  n-n  , 


=  0  otherwise. 
We  then  obtain 


x  ■  5  + 


t  I  +  ?  I  "J 


(26) 


n-n 


where  c  =  no.  of  nearest  neighbors  per  atom  , 

A  -  & 

A  "  2  ' 


Substituting  (23)  and  (26)  into  (25)  yields 

/E  +|i 


=  ex 


=  exp 


<{  (4  ’1)41  •»  E  v,  *  <■.«*$  *  E  *.} 


(o^) 


n-n 


{(4  4)  si} 


(27) 


where  Z  is  the  Ising  partition  function  of  eq.  (2)  provided  we  make  the 
following  identifications: 


J  = 


H  =  (2E  +  2\x  +  cV)  /  h 

s 


(28a) 

(28b) 


It  is  also  clear  that  the  fractional  converage  of  the  surface  6  is 
related  to  the  magnetization  per  spin  M  by 


■ 


. 
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M  -  2  0-1 


(28c) 


When  the  absorbed  phase  is  in  equilibrium  with  a  gas  phase  the 
chemical  potentials,  |i,  of  the  two  phases  are  equal.  For  a  classical 
ideal  gas  the  chemical  potential  is  given  by  (see,  for  example,  ref.  39) 


e“>T  ,  (  h*  V/g 

V  2tt  m  kT  ) 

where  h  is  Plank's  constant, 
p  is  the  gas  pressure, 
and  m  is  the  mass  of  a  molecule. 


(29) 


The  case  of  H  =  0  corresponds  to  a  particular  value  of  pressure 
at  each  temperature.  It  is  convenient  to  introduce  certain  parameters. 

E 


Let 


€  = 


t  = 


a 


V 

T_ 

Tc 

J 

kTc 


(interaction  ratio)  , 

(dimensionless  temperature), 

y “TT"  (  =  \  sinh  .  1+U07  for  a  square 

4  d  lattice)  , 


=  (2nm)^2  V^2  (ha)  h  ^  (base  pressure)  , 


P  =  (dimensionless  pressure) 


With  these  definitions  (29)  becomes 
e^T  .  P  t'5/2  . 


The  condition  for  H  =  0  becomes 

P  =  Po  =  t5/2  e-(U+2c)a/t  _  (30) 

At  this  pressure  we  will  have,  through  the  correspondence  with  the  Ising 
ferromagnet,  0  =  ^  if  t  >  1  and  a  discontinuity  in  0  corresponding 


' 
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to  spontaneous  magnetization  if  t  <  1 


The  coverage  0  will  be  a  function  of  the  two  variables  (of  eq.(3)) 
-2a 


and 


x  =  e 


y  =  P" 


(30 


The  symmetry  of  the  Ising  p.f.  under  magnetic  field  reversal, 
M(H)  =  -  M(-H)  ,  induces  a  symmetry  to  the  adsorption  isotherms 


s(  p-,  0  -  i  -  e(  p2,  t) 


(32) 


A  final  correspondence  which  is  useful  is  between  the  ferromagnetic 
susceptibility  and  the  slope  of  the  isotherm.  At  H  =  0  we  have 


39 

d(f-) 

o 


P  =  P 


3m 

3v 


=  r  x  («) 


(33) 


V  =  0 


where  x(u)  Is  called  the  reduced  susceptibility  at  zero  field  and  is 
related  to  the  physical  susceptibility  by 


3m 

5h 


H  =  0  "  kT 


Using  these  relations  it  is  possible  to  calculate  very  accurate 
adsorption  isotherms  for  this  model. 


9-  Series  Expansions  for  the  Plane  Square  Ising  Lattice. 

The  method  of  calculating  the  expansions  for  In  A  and  in  ft  is 
based  on  eq.  (10),  (11)  and  (12).  It  illustrates  the  use  of  symmetries  and 
special  properties  which  may  be  known  in  evaluating  the  p.f,. 


- 


. 


- 
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The  form  of  the  calculation  is  given  in  table  II.  We  begin  with 

the  values  of  the  coefficients  g(N,Y)  in  in  A  for  V  <  1U  which  have 

been  found  by  Domb  [30].  Only  a  finite  number  of  these  are  non-zero  since 

in  a  lattice  which  is  infinite  in  more  than  one  direction,  larger  values 

of  N  ultimately  require  larger  values  of  Y.  The  coefficients  for 
been 

Y  <  12  have^counted  (Domb  used  a  matrix  method)  as  a  check. 

In  any  lattice  for  which  the  coordination  number  c  is  an  even 
constant  it  is  easy  to  show  that  g(N,Y)  =  0  for  Y  odd.  The  proof  by 
induction  fallows  by  observing  that  overturning  an  additional  spin  will 
create  c  (if  isolated),  c  -  2  (if  adjacent  to  one  other  (})  spin), 
c  -  k,  etc.  (if)  pairs.  Thus  the  parity  of  Y  can  not  change.  The 
statement  is  clearly  true  for  N  =  1.  If  c  is  an  odd  constant  then  the 
same  argument  shows  that  Y  always  changes  parity  when  N  is  increased 
by  one.  So  we  would  then  have  g(N,Y)  =  0  unless  N  and  Y  are  both 
even  or  both  odd. 

In  order  to  find  the  coefficients  i,m )  we  write  eq.  (10) 
with  them  as  unknowns,  making  use  of  the  known  values  of  g(N,Y).  Init~ 
ially,  we  have  more  equations  than  unknowns.  The  existence  of  a  solution 
to  such  a  set  provides  a  check  on  the  values  of  the  coefficients  with  which 
we  started.  When  the  number  of  unknowns  reaches  the  number  of  equations 
we  reverse  the  procedure  and  regard  the  g(N,Y)  as  unknowns,  using  the 
newly  found  values  of  2rlt  ,m).  However  the  procedure  does  not  go  on 
indefinitely,  as  the  number  of  unknowns  increases  more  rapidly  than  the 
number  of  available  equations.  We  see  that  i  =  N  =  7  is  as  far  as  we 
may  procede.  The  values  of  are  now  calculated  using  eq.  (ll). 

This  method  is  very  flexible,  however,  and  any  available  inform¬ 
ation  about  the  coefficients  may  be  used  to  extend  it.  For  example  we  may 
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determine  Jb  ( i ,0)  independently  by  means  of  the  duality  theorem  [I13], 


(31*) 


N 

This  enables  us  to  check  our  previous  values  of  £  (i,o)  and  to  extend  the 
calculations,  using  eq.  (ll),  to  obtain  (8, m) ,  g(8,Y),  and  J?/(9,m)  . 

Sykes  [3U]  has  found  several  values  of  J^(i,2)  by  means  of 
his  counting  theorem.  Using  his  values  the  calculation  may  be  extended 
to  obtain  ( 10, m) ,  g(9»Y)»  and  ^(ll,m). 

It  is  clear  from  this  example  that  the  exploitation  of  symmetries 
and  special  properties,  using  (10),  (ll),  and  (12),  can  be  a  very  useful 
and  powerful  technique. 

The  calculation  described  above  has  only  been  carried  out  up  to 
i  =  N  =  8,  since  this  is  sufficiently  accurate  for  our  present  application 
and  the  labor  involved  increases  with  the  order  of  the  coefficients.  The 
values  of  g(N,Y)  obtained  agree  with  those  of  Domb  and  Sykes  [32]. 

(Domb  and  Sykes  give  the  series  for  A  and  we  must,  rather  laboriously, 
take  its  logarithm  in  order  to  compare.  Historically  this  was  done  first 
and  it  was  later  verified  that  the  above  equations  were  satisfied  when  the 
’'unknowns"  were  g(N,Y).)  Ref.  [18]  contains  a  table  of  graphs,  classified 
topologically,  from  which  we  may  find  Jb  (i,m)  for  i  <  5  •  Our  calculat¬ 
ed  values  are  in  agreement  with  these. 


The  series  for  in  A  and  in  Q  which  we  have  calculated  are 


listed  in  appendix  III. 
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TABLE  II  -  FORM  OF  THE  CALCULATION  FOR  A  ONE  LAYER  SQUARE  LATTICE 

c  =  4,  7=2 

g(N,Y)  =  0  for  Y  odd. 

Given  g(N,Y)  for  Y  <  14. 


Unknowns 

No. 

Equations 

No. 

Extra 

Equations 

1 

(10)  1=1,  N-1...4 

4 

3 

^(2,m) 

2 

(10)  1=2,  N=1...4 

4 

2 

3 

(~10)  4=3,  N=1...4 

4 

1 

i/(4,m) 

4 

(10)  4=4,  N=1...4 

4 

- 

8(5, V) 

3 

flO)  N=5,  4=1.'.".  4;  (l2) 

5 

2 

5 

(10)  4= 5,  N=l. . .5 

5 

- 

g(6.Y) 

5 

(10)  N=5,  1=1... 5;  (12) 

6 

1 

6 

(10)  4=6,  N=1 ...6 

6 

- 

g(7,V) 

7 

(10)  n=7,  1=1... 6;  (12) 

7 

- 

7,tn) 

7 

(10)  4=7,  N=l...7 

7 

- 

i^(8,m) 

8 

(10)  4=8,  N=l. . .7;  (11)* 

8 

- 

s(8,Y) 

9 

(10)  N=8,  4=1... 8;  (12) 

9 

- 

9 

(10)  4=9,  N=l...8;  (11) 

9 

- 

^f(lO,m) 

9  r 

(10)  i=io,  n=i...8;  (li)* 

9 

- 

g(9,Y) 

11 

(10)  N=9,  4=1... 10;  (12) 

11 

- 

i? (11, m) 

10  t 

(10)  4=11,  N=l. . .9;  (11) 

10 

- 

# 


t 


Using  the  duality  theorem  (^4)  • 
Using  Sykes'  [34]  value  of  2)  . 
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10.  Calculation  of  the  Isotherms, 

The  qualitative  behavior  of  the  isotherms  is  a  smooth  increase 
of  0  with  increasing  P  for  t  >  1.  However  for  t  <  1  there  is  a 
finite  jump  in  0  at  a  particular  pressure  P  =  ,  which  increases 

with  t.  This  phenomenon  is  known  as  critical  adsorption. 

Because  P  =  P^  corresponds  to  H  =  0,  it  is  possible  to 
calculate  the  critical  region  of  the  isotherms  exactly  using  Yang’s  [26] 
formula  for  the  spontaneous  magnetization, 

M(0)  =  [ ( 1+x2)  (l-x2)~2  ( l-6x2+x^)^  ]4  (35) 

where  x  =  e 

This  is  interesting  since  approximate  calculations  are  frequently  un¬ 
reliable  in  the  critical  region.  From  equations  (28c),  (30),  and  (35) 
we  have  plotted  in  Fig.  1  the  curves  of  coverage  0q  versus  transition 
pressure  P^  for  two  values  of  interaction  ratio  e  .  The  main  effect 
of  changing  c  is  to  change  the  pressure  scale,  but  a  change  of  shape 
is  also  displayed. 

To  calculate  the  isotherms  for  P  /  we  use  equations  (18) 

and  (22)  with  (28),  (31)  and  the  series  expansions  found  in  section  9- 

It  is  known  from  the  work  of  Lee  and  Yang  [ UU ]  that  the  series  for  in  A 

converges  for  j y |  <  1  if  0  <  x  <  1  .  Thus  we  obtain  6  for  all  t 

and  for  P  >  P  .  By  the  symmetry  property  (32)  we  also  have  0  for 
o 

P  <  P  .  Since  the  exact  values  of  0  at  P  =  P  are  known,  we  have 
o  o 

here  a  fairly  good  calculation  of  the  complete  isotherms. 

For  t  >  1  the  in  0,  series  is  much  better  than  in  A  near 

P  ^  P  .  Less  is  known  about  the  convergence  conditions  of  the  in  Q 
o 

series.  It  converges  for  ju|  <  uc  =  \T~ 2  -  1  when  v  =  0  , 


2J/kT 


and  on 


tw  ■ 
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the  basis  of  numerical  calculations  the  convergence  radius  seems  to  be  at 
least  this  large  for  -  1  <  v  <  1  .  This  then  gives  us  isotherms  for 
t  >  1  and  all  pressures. 

Using  (33)  we  can  obtain  the  slope  of  the  isotherm  at  P  =  P 

o 

for  t  >  1  from  Sykes>  [3^]  expansion  of  the  ferromagnetic  susceptibility, 

15 

which  he  has  found  up  to  u  .  Calculating  the  slope  from  our  In  to 
series  would  be  equivalent  to  using  the  first  eight  terms  of  Sykes’  series, 
so  this  provides  a  check  on  the  round  off  error  committed. 

A  further  check  on  the  calculations  is  provided  by  the  fact  that 
the  two  series.  In  A  and  In  to  ,  both  converge  in  the  range  t  >  1, 

0  <  y  <  1  and  can  be  compared  in  this  region. 

The  results  of  these  calculations  are  displayed  in  figure  2. 

By  examining  the  behavior  of  the  coefficients  in  the  series,  comparing 
the  sums  of  the  two  series  with  each  other  and  with  the  exact  values  at 
P  =  ,  we  estimate  the  accuracy  of  this  treatment  as  follows.  The 

error  in  9  is  less  than  0.01  for  the  ranges: 

t  <  0.9  or  t  >  1.6  and  all  P  ; 

§  <0.7  or  p  >  1.^  and  t  . 
o  o 

For  the  remaining  region  the  error  may  be  somewhat  larger,  but  it  can  not 
be  very  great  since  at  P  =  PQ  ,  ©  is  known  exactly  and  the  slope  of 

the  isotherm  is  known  to  a  high  accuracy  for  t  >  1  . 
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Fig.  1  Coverage  at  transition,  c  , versus  dimensionless 
transition  pressure,  . 
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Pig.  2  Theoretical  Adsorption  Isotherms 


-  27  - 


11 .  Comparison  with  Experiment* 

Because  of  the  simplifying  assumptions  of  our  model  such  as  a 
regular  array  of  adsorption  sites,  monolayer  coverage  only,  and  no  internal 
degrees  of  freedom,  a  comparison  of  the  theoretical  isotherms  of  fig.  2 
with  experimental  isotherms  should  not  be  taken  too  seriously.  However,  as 
mentioned  in  Section  2,  the  data  of  ref.  [6],  [7],  and  [8]  exhibit 
transitions  very  similar  to  our  "ferromagnetic"  monolayer  transitions. 

We  shall  therefore  procede  to  analyse  their  data  in  terms  of  our  model. 

A  knowledge  of  the  critical  temperature,  T  ,  and  critical 
pressure,  p^,  above  which  the  sudden  transitions  in  9  no  longer  occur 
is  sufficient  to  determine  the  adsorption  isotherms  for  our  model.  First 
V  follows  immediately  from  .  Then  V  and  m  determine  the  base 

pressure  p6.  From  the  critical  pressure  one  gets  Pc  =  Pcy/P°  •  Finally 
from  (JO)  we  obtain 

€  =  -  2  -  0.56728  in  Pc 

and  E  =  €  V  . 

s 

The  information  thus  derived  from  the  experimental  curves  of 
ref.  [63s  [7],  and  [83  is  tabulated  in  table  III.  The  values  of  V 
and  E  are  the  right  order  of  magnitude.  Also  Eg  depends  mainly  upon 
the  substrate,  as  one  would  expect  physically. 

• 

It  is  next  possible  to  compute  complete  sets  of  isotherms  for  the 
six  cases.  We  have  chosen  to  do  this  only  for  krypton  on  sodium  bromide 
because  in  the  case  of  n  =>  heptane  on  ferric  oxide  the  coverage  appreciably 


*  The  calculations  in  this  section  were  done  by  Dr.  Betts. 
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excedes  a  monolayer  for  pressures  only  slightly  above  critical  and  in  the 
other  cases  the  data  is  less  complete. 

Theoretical  and  experimental  isotherms  are  compared  in  fig.  3. 
By  using  experimental  values  of  and  p^  we  have  essentially  fitted 

the  theoretical  curves  at  the  point  T  =  77.02°  K,  p  =  0.112  mm.  , 

6  =  0.5*  Also  rather  than  accept  Fisher  and  McMillan's  fractional  cover¬ 
ages  we  have  miltiplied  their  0  values  by  0.875  so  that  the  midpoint 
of  the  discontinuity  is  at  Q  =  0.5*  as  required  by  our  model. 

The  resulting  sets  of  isotherms  are  qualitatively  similar. 
However,  the  experimental  variation  of  transition  pressure,  p^,  with 
temperature  is  less  than  the  theoretical  and  the  experimental  magnitude 
of  the  discontinuity  in  6  is  less  than  the  theoretical. 


The  former  discrepancy  is  probably  due  to  the  neglect  of 
vibrational  degrees  of  freedom  of  the  adsorbed  molecules.  It  can  be 
shown  thermodynamically  that  (see  appendix 

„  «  t(CPS-V/R  -L./H 

p  c£  1  e  ' 

o 


where  C  and  C  ^  are  the  specific  heats  (  at  const,  pressure)  of 

the  gaseous  and  adsorbed  phases  respectively s  and  Lo  is  the  latent 

heat  of  the  phase  transition  at  T  =  0  .  Putting  C  -  ~  and  C  -  0 

Pg  2  pa 

we  recover  essentially  eq.  (30).  Thus  by  postulating  appropriate  surface 
degrees  of  freedom  we  could  formally  explain  the  former  discrepancy. 


The  second  discrepancy  may  be  due  to  imperfections  in  the  sub¬ 
strate  lattice.  The  effect  of  these  would  be  to  hinder  the  propagation  of 
long  range  order,  which  would  result  in  a  less  sudden  setting  in  of  spont¬ 
aneous  magnetization  in  the  corresponding  model  of  a  ferromagnet.  Thus 
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the  magnitude  of  the  discontinuity  at  a  given  t  would  be  smaller. 

We  can  be  quite  sure  that  the  type  of  lattice  chosen  (square) 
is  not  responsible  for  the  latter  discrepancy  since  the  spontaneous  mag¬ 
netization  of  the  triangular  and  hexagonal  lattices  differ  very  little 
from  that  of  the  square  lattice  (see  ref.  [18],  p.  217,  Fig.  32). 

The  fact  that  adsorption  is  not  really  restricted  to  one  layer 
may  also  play  some  part  in  the  latter  discrepancy.  It  is  difficult  to  say 
just  what  its  effect  would  be  in  the  critical  region,  except  that  the 
symmetry  property  (32)  would  cease  to  be  valid. 

We  conclude  by  emphasizing  that  any  substantial  disagreement 
between  theory  and  experiment  here  is  due  to  inadequacies  of  the  model 
and  not  to  any  mathematical  approximation. 
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IV  THE  MANY  LAYER  ISING  MODEL 

12.  Introduction. 

To  avoid  confusion  we  should  point  out  that  the  n- layer  Ising 
model  does  not  provide  the  generalization  of  the  treatment  in  the  previous 
chapter  to  the  case  of  multilayer  adsorption.  This  is  so  because  the 
summation  over  all  cr^  =  +  1  would  include  configurations  in  which 
molecules  in  higher  layers  were  "adsorbed"  on  top  of  holes  in  the  first 
layer.  However  this  model  is  of  interest  for  two  reasons.  First,  as 
already  stated  in  Section  2,  it  is  a  simple  model  of  a  cooperative  phen¬ 
omenon  in  a  thin  film  which  exhibits  a  shift  of  A  -  point  with  the  number 
of  layers.  This  is  qualitatively  similar  to  a  number  of  experimental 
results,  [11,12,13]. 

Secondly,  Domb  throughout  his  review  article  [18]  expresses  the 
view  that  certain  features  of  the  Ising  model  are  dependent  only  or  prim¬ 
arily  upon  the  dimensionality  of  the  lattice.  Since  the  n-layer  lattice 
is,  in  some  sense,  intermediate  between  two  and  three  dimensions,  it  is  of 
theoretical  interest  to  see  how  the  properties  vary  with  n  . 

It  was  hoped  that  the  n-layer  lattice  might  also  be  intermediate 
in  difficulty  between  the  two  and  three  dimensional  lattices,  and  that  an 
accurate  determination  of  the  variation  of  critical  temperature  with  n 
would  be  possible.  Unfortunately  this  proved  not  to  be  so  (with  the 
exception  of  n  =  2)  and  only  a  crude  determination  was  possible. 

In  this  chapter  we  shall  mean  by  a  "cubic  lattice"  any  of  the 
n  -  layer  lattices  whose  primitive  unit  cell  is  a  cube.  The  plane  square 
lattice  can  be  considered  a  one  layer  cubic.  The  usual  simple  cubic 
lattice  will  be  called  the  infinite  cubic. 


;  ' 
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13 •  Calculation  of  the  Partition  Functions. 

The  calculations  of  the  series  expansions  for  2,  3,  and  4  layers 
are  similar  to  that  for  one  layer,  which  was  explained  in  Section  9.  The 
tables  IV  and  V,  which  outline  these  calculations,  are  self  explanatory. 

A  number  which  measures  the  degree  of  difficulty  which  a  parti¬ 
cular  lattice  presents  for  this  method  is  the  rate  at  which  unknowns  in¬ 
crease  relative  to  the  rate  at  which  equations  become  available.  Spec¬ 
ifically,  consider  the  unknowns  in  eq.  (lo)  to  be  g(N,Y).  Then  as  N  is 
increased  by  one  the  number  of  unknowns  increases  by  c.  For  a  homogeneous 
lattice  (one  in  which  all  sites  are  equivalent)  we  showed  in  Section  9  that 
g(N,Y)  =  0  for  Y  odd  if  c  is  even,  and  g(N,Y)  =  0  for  N  and  Y  of 
opposite  parity  if  c  is  odd.  Thus  half  of  the  unknowns  are  zero  for  a 
homogeneous  lattice. 

We  can  therefore  say  the  "difficulty"  of  the  square  lattice  is 

c 

—  =  2  .  The  two  layer  cubic  is  a  homogeneous  lattice  with  c  =  5  •  Its 
"difficulty"  is  .  The  "difficulty"  of  the  infinite  cubic  is  3*  How¬ 
ever  the  n-layer  cubic  for  2  <  n  <  00  is  inhomogeneous,  since  surface  sites 
are  not  equivalent  to  interior  sites.  Therefore  its  "difficulty"  is  c  =  6. 
The  greater  the  "difficulty"  the  lower  the  order  to  which  coefficients 
may  be  calculated  by  our  method. 

It  is  interesting  to  note  that  the  linear  chain  has  "difficulty" 
one,  i.e.  the  unknowns  increase  at  the  same  rate  as  the  equations.  This 
means  that  the  coefficients  up  to  any  order  can  be  found  and  so  the  part¬ 
ition  function  is  completely  determined  by  eq.  (10),  (ll)»  and  (12). 

Although  the  n-layer  cubic  proved  to  be  more  difficulty  than  the 
infinite  cubic  in  terms  of  applying  equations  (10),  ( 1 1 ) ,  and  (12),  we  have 


' 


■  1 

. 
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TABLE  IV  -  FORM  OF  THE  CALCULATION  FOR  TWO  LAYERS 

c  *=  5*  7  =  5/2 

We  have  counted:  g(N,Y)  for  Y  <  I7  , 

J?/(4,0)  for  4  <  10  . 

g(N,Y)  =  0  unless  N  and  Y  are  both 
odd  or  both  even  . 


Unknowns 

No. 

Equations 

No. 

Extra 

Equations 

Jb  (l^tt) 

1 

(10)  4=1,  N= 1 , 2 , 3 

(11) 

4 

3 

&  (2,m) 

2 

(10)  4=2j  N=  1 , 2 , 3 

(n) 

4 

2 

3 

(JO)  <=5,  N=l,2,3 

(11) 

4 

1 

4 

(10)  i»4,  N= 1 , 2 , 3 

(n) 

4 

- 

g(4,Y) 

2 

(10)  N=4,  1=1... 4 

(12) 

5  . 

.  3 . 

g(5,V) 

4 

(10)  N=5,  i=l. . .4 

(12) 

5 

1 

^(5,^) 

5 

(10)  4=5,  N=1...5 

(u) 

6 

1 

jy(6,m) 

6 

(10)  4=6,  N-1...5 

(u) 

6 

- 

g(6,Y) 

7 

(10)  N=6,  4=1... 6 

(12) 

7 

- 

J*(7,m) 

7 

(10)  4=7,  N-1...6 

(11) 

7 

- 

TABLE  V  -  FORM  OF  THE  CALCULATIONS  FOR  5  AND  4  LAYERS 
c  =  6,  7=3“  l/n  ,  n  =  no.  of  layers 

We  have  counted:  g(N,Y)  for  Y  <  12  and  all  N, 

N  <  3  and  all  Y, 

0)  for  4  <  10. 


Unknowns 

No. 

Equations 

No. 

Extra 

Equations 

J*(  l»m) 

1 

(10)  4=1,  N-1,2,3; 

(11) 

4 

5 

J??(2,m) 

2 

(10)  4=2,  N=1,2,3J 

(iiT 

4 

2 

M  3,®) 

3 

(10)  4=3,  H-1,2,3; 

(11) 

4 

1 

-M4,m) 

4 

(10)  4=4,  N=l,2,3; 

(U) 

4 

- 
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developed  a  method  by  which  closed  graphs  in  any  n-layer  cubic  may  be 
easily  found  after  they  have  been  enumerated  by  space  type*  for  the 
infinite  cubic.  We  observe  that  the  number  of  ways  in  which  a  particular 
space  type  can  be  oriented  in  n  layers  is  smaller  than  in  the  infinite 
cubic  because  of  the  finite  distance  to  the  surface.  We  can,  as  we  shall 
show,  relate  the  contribution  of  a  space  type  in  n  layers  to  its  contri¬ 
bution  in  the  infinite  cubic  by  a  relation: 

i/n  "  fn  (a.b.c)  Jt> „  (57) 

on 

where  fn(a,b,c)  is  called  an  orientation  factor  and  depends^the  number 
of  layers  n  and  the  dimensions  (a,b,c)  of  the  smallest  rectangular 
box  which  will  enclose  a  graph  of  the  space  type. 


Now  the  number  of  ways  in  which  a  box  of  dimensions  (a,b,c)  can 
be  oriented  in  three  dimensional  space  by  rotations  of  the  cubic  group  is 
24.  If  it  is  inside  an  n=layer  cubic  the  number  of  orientations  per  site 
with  the  a-axis  perpendicular  to  the  plane  of  the  layers  will  be 


8 


(n~a) 


n 


for  n  >  a  ,  zero  for  n  <  a  , 


and  similarily  for  the  b  and  c  axes. 


Thus 


fnU,b,c) 


1_ 

2b 


Sinral  8_{,n:b),  +  8(.r tlgl 

n  n  n 

s=  (3n-a»b“c)  /  3n  if  n  >  a,  b,  c 
f  (a,b,c)  =  (2n-b"c)  /  3n  if  a  >  n  >  b,  c  , 
fn(a,b,c)  =  (n~c)  /  3n  if  a,  b  >  n  >  c  , 


(38) 


fn(a,b,c) 


=  0 


if  a,  b,  c  >  n 


Graphs  equivalent  under  the  48  operations  of  the  cubic  group  form  a 
space  type.  The  term  is  essentially  synonomous  with  shape. 


J 
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By  means  of  equations  (37)  ant^  (3$)  we  may  now  find  the  number  of  connected 
graphs  which  can  be  drawn  on  an  n-layer  lattice. 


For  disconnected  graphs  we  must  use  a  method  due  to  Domb  [18]. 

Allow  the  component  figures  to  range  independently  over  the  lattice  and 
classify  the  graphs  which  will  result.  For  two  squares: 

(  D)2  =  2(D  □  )  +  2(00  )  +  s(oa  )  +  (a  )  , 

where  the  pictorial  symbols  indicate  the  number  of  graphs  of  that  topology 
which  can  be  drawn  on  the  lattice.  Thus  the  disconnected  graphs  are  express¬ 
ible  in  terms  of  connected  graphs  to  which  the  above  method  of  orientation 
factors  is  applicable.  Formulas  for  disconnected  configurations  are  given 
in  ref.  [18],  p.  356. 


Those  used  in  this  thesis  are 

(on)  +  (00)=  §(  a  )s  -  (m  )  -  |(  o )  (uoa) 


(o  □)  +  (OO)  =  (o)  x  (a)  -3(ag)  )-(cO  )-2«2>  )-2(m  )  (4ob) 


=  B2  P,  P 


6  b 


b(3P 


LOb 


9k 


2P 


8c 


2P  ) 


The  literal  notation  is  Domb ' s  and  corresponds  to  the  pictorial  notation 
above.  Domb's  formulas  also  include  some  topological  types  which  can  not 
occur  on  a  cubic  lattice  and  so  have  been  omitted  above. 


In  table  VI  we  have  classified  by  dimensions  all  the  graphs 
necessary  to  evaluate  ,0;B)  for  £  <  10  on  an  n-layer  cubic.  Our 

values  of  Jr«  ,0;B)  for  the  infinite  cubic  agree  with  those  of  Rushbrooke 


and  Eve  [35]. 


. 
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This  method  of  orientation  factors  is  applicable  to  any  type  of 


graph  but  it  becomes  impractical  for  other  than  closed  graphs  because  of 
the  much  larger  number  of  space  types  which  must  be  enumerated.  For 
example,  there  are  3  space  types  of  hexagons  and  11  space  types  of  octagons, 
but  there  are  at  least  I77  and  4042  space  types  of  loose  ended  chains  with 
six  and  eight  bonds  respectively. 


The  reduced  susceptibility  at  zero  field  is  found  from  (22)  and 

(33)  to  be 


00 


(41) 


1=1 


Studying  the  terms  of  this  expansion,  for  homogeneous  lattices,  Sykes  [34] 
has  found  a  powerful  "counting  theorem",  the  only  serious  flaw  of  the 
"theorem"  being  that  it  lacks  a  general  proof.  However  in  the  process  of 
establishing  the  "theorem"  up  to  order  eight,  Sykes  has  derived  some  useful 


formulas.  We  denote  i/(i,2)  =  , 


and 


s 


Then 


(U2) 


A  •  2sA-i  +  s2  A- 2  =  Al  '  1  *  5  ’ 


where  for  a  cubic  lattice 


■ 
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d 

d 

d 

d 

d 

d 

d 


3 

4 

5 

6 

7 

8 

9 


0 


-4  P, 


0 


4PU  -  6P6 


(43) 


b  P. 


'7® 


4  %  +  6  P6  +  16  V  '  8  p8  +  4  P8c  +  4  P8h 


11  V  +  8  P8c  +  8  P8h  +  4  P9j  +  4  V  +  4  % 


The  last  term  is  given  by  Domb  (ref.  [18],  p.  324). 
equal  4  of  Domb's  and  Sykes'  d  . 

C.  Jb 


Note  that  our  d 

Ju 


Using  (43)  we  have  derived  ^(8,2)  and  >^/(9,2)  for  the  two 
layer  lattice.  The  graphs  needed  in  (43)  are  classified  in  table  VI.  We 
prefer  not  to  use  the  most  powerful  form  of  the  counting  theorem  since  it 
has  not  been  proven.  Although  we  did  in  chapter  III  evaluate  numerically 
the  fifteen  term  susceptibility  series  for  the  square  lattice,  which  was 
derived  from  the  counting  theorem,  it  has  been  independently  derived  up 
to  the  eleventh  term.  So  if  the  theorem  should  be  false  we  would  only 
have  to  increase  our  error  estimate  slightly. 


Since  eq.  (42)  explicitly  involves  the  coordination  number  c, 
it  is  difficult  to  see  how  these  results  may  be  generalized  to  an  inhomo¬ 
geneous  lattice.  Thus  here  again  the  n- layer  cubic  (2  <  n  <  00)  appears 
more  difficult  than  the  infinite  cubic. 


The  series  which  have  been  derived  in  this  section  are  listed 


in  appendix  III. 


m  ■ 


■ 
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TABLE  VI  -  GRAPHS  ON  A  SIMPLE  CUBIC  LATTICE 


Simple  Polygons 

Composite  Figures 

dimensions 

contribu- 

dimensions 

contribu- 

a,b  ,c 

tion  to 

a,b ,  c 

tion  to  Goo 

P4 

1,1,0 

3 

V  > 

2,1,0 

6 

P6 

2,1,0 

6 

1,1,1 

12 

1,1,1 

16 

P8c  (<2>  > 

1,1,1 

24 

P8 

3,U0 

6 

p8h  > 

2,2,0 

6 

2,2,0 

15 

2,1,1 

24 

2,1,1 

180 

p9k  (C0  > 

3,1,0 

12 

1,1,1 

6 

2,2,0 

24 

P10 

4,1,0 

6 

2,1,1 

264 

3,2,0 

78 

1,1,1 

24 

3,1,1 

408 

P9j  (O-O) 

3,2,0 

12 

2,2,1 

1668 

3,1,0 

12 

2,1,1 

252 

3,1,1 

24 

2,2,1 

60 

2,1,1 

72 

1,1,1 

12 

V10 ) 

2,1,1 

12 

1,1,1 

8 

W08  > 

2,1,1 

12 
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111.  Results . 

Using  the  series  derived  in  the  previous  section  with  equations 
(I7)  and  (21),  we  have  plotted  in  figure  1|  the  specific  heats  at  zero  field 
for  2,  3,  and  4  layer  Ising  lattices.  The  curves  for  n  =  1  and  n  =  00 
come  from  Onsager's  exact  solution  for  the  square  lattice  [20]  and  Wake¬ 
fields  series  expansions  for  the  simple  cubic  [31].  We  see  that  critical 
temperature  seems  to  vary  smoothly  with  n  . 


In  order  to  estimate  the  critical  temperature  more  accurately  we 
apply  a  method  due  to  Wakefield  [45].  The  critical  temperature  is  marked 
by  a  singularity  in  the  p.f..  This  means  that  a  power  series  in  some 
temperature  variable  (eg.  u  =  tanh  J/kT)  will  have  a  finite  radius  of 
convergence.  We  therefore  wish  to  estimate  the  radius  of  convergence  of 
certain  power  series. 


Consider  a  series 


^  a^  u11  ,  convergent  for  |u|  <  u^  .  By 


the  ratio  test 


n=0 


lim 
n  — *■<» 


a 

n 


1 

u 

c 


Then 


Let  us  assume 


P 

a  A  —  for  large  n  . 

n  n 

u 

c 


(1  ♦  Y 

n 


(MO 


1  d  +  £) 

u  n' 

c 


(45) 


■  ■  ■; 


37  - 


Thus  u  and  P  can  be  estimated  from  a  plot  of  -  versus  —  .  The 

c  a  ,  n 

n-1 

index  P  is  related  to  the  nature  of  the  singularity  as  follows  [46] I  for 
P  not  a  negative  integer,  the  singularity  is  of  the  form  (l  -  Z)"(P+1)  > 
z  »  u/iic  ;  for  P  a  negative  integer,  the  singularity  is  of  the  form 
( 1-z)  ^  ^  /n(l-z)  .  Thus  the  position  and  nature  of  a  singularity  can  be 

estimated  from  the  coefficients  of  a  power  series. 


There  are  two  series  to  which  we  may  apply  this  method.  If  we 
choose  a^  =  J??(/,0)  we  shall  be  considering  the  zero  field  p.f.,  and  will 
be  looking  for  the  singularity  in  internal  energy  and  specific  heat.  If  we 
choose  a^  =  h(t,2)  we  shall  be  considering  the  zero  field  susceptibility. 
The  former  approach  has  been  used  by  Wakefield  [45]  and  the  second  by  Domb 
and  Sykes  [36].  Although  no  mathematical  proof  is  known,  one  can  argue 
physically  that  the  two  series  should  yield  the  same  value  of  uc  .  The 
reason  is  that  the  setting  in  of  long  range  order  gives  rise  to  the  sing¬ 
ularity  in  internal  energy  and  to  spontaneous  magnetization,  which  implies 
the  susceptibility  becoming  infinite. 


For  the  p.f.  series  our  data  gives  only  three  ratios,  which  hardly 
justifies  any  extrapolation.  The  susceptibility  series  is  no  better  for 
n  =  3  or  4,  but  we  have  several  terms  for  n  =  2.  The  ratios  are  plotted 
in  figure  5.  Our  estimate  of  the  parameters  for  n  =  2  are  given  in  table 
VII  below,  along  with  the  Domb's  and  Sykes*  values  for  n  =  1  and  00  . 

Table  VII  -  Asymptotic  Parameters  for  Susceptibility  Series 


n  u  P 

c 


1 

0.4142 

0.75 

2 

0.301 

0.74 

00 

0.2182 

0.25 

approximate  variation  of 

Xc  =  e"2J/kTc  with  n  is  shown 

in  fig.  5a. 


n  is  shown 
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Figure  5  Upper  curves  »  lower  curves 

Jo  \  "  *  s>  ) 


—  [  versus  7  for  2  0  ,  3  A  ,  and  4  H 

L  ^(i-2,0)  -I  i 


layer  lattices, 


.  2J  /kTc 

Figure  5a  =  e  x  versus  number  of  layers,  n 

The  vertical  bar  for  n  =  3  an^  4  emphasizes  the 


uncertainty  but  does  not  represent  confidence  limits. 
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V  MATHEMATICAL  QUESTIONS  RAISED 
The  work  of  this  thesis  has  suggested  certain  mathematical 
questions.  The  first  and  probably  the  simplest  concerns  equation  (H.T.,E) 
of  table  I.  This  equation,  giving  the  total  number  of  graphs  with  a  given 
number  of  odd  vertices  which  can  be  drawn  on  a  lattice,  should  be  derivable 
from  the  properties  of  graphs  without  making  use  of  the  Ising  partition 
function . 


Secondly,  equation  (10)  relates  the  coefficients  g(N,Y)  and 
M*  ,m).  But  these  two  sets  of  coefficients  are  both  determined  by  the 
topology  of  the  lattice.  It  should  therefore  be  possible  to  prove  (10) 
from  topological  considerations  without  the  use  of  the  Ising  p.f..  Such 
a  proof  would  give  greater  insight  into  the  meaning  of  eq.  (lo)  and  might 
provide  a  starting  point  for  further  progress  in  the  theory  of  cooperative 
phenomena.  It  is  significant  to  note  that  Sykes'  unproven  counting  theorem 
was  derived  from  a  topological  analysis  of  the  graphs  involved. 

The  third  concerns  the  consequences  of  magnetic  field  reversal 
invariance  on  the  coefficients  of  the  L.T.  series  for  an  infinite  lattice. 
The  transformation  H  — ►  -H  corresponds  to  v  -*•  -v  in  the  H.T.  series. 

Thus  invariance  simply  means  the  absence  of  odd  powers  of  v,  i.e. 

JfU  ,m)  =  0  for  m  odd.  The  corresponding  transformation  for  the  L.T. 
series  is  y  — ►  l/y  .  From  the  work  of  Lee  and  Yang  [kU]  we  know  that 
all  the  singularities  of  in  A  (x»y)  in  the  complex  y-plane  lie  on  the 
unit  circle  (0  <  x  <  l)>  and  above  the  critical  temperature  they  do  not 
form  a  dense  boundary  in  the  neighborhood  of  the  real  axis.  Thus  analytic 
continuation  beyond  the  unit  circle  is  possible.  So  we  ask  the  following 
question,  which  is  a  slight  simplification  of  question  appropriate  for 


the  p.f.. 


Hr* 1 '  *v 

- 


aim  mm 
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Consider  the  power  series 
00 

f(y)  -  I  an  y»  ,  convergent  for  |y|  <  1  . 
n=0 

What  are  the  necessary  and  sufficient  conditions  on  the  coefficients  a 

n 

such  that  the  analytic  continuation  of  f(y)  satisfies  the  identity 

f(y)  =  f(^)  ? 

For  the  Ising  p.f.  we  identify  f(y)  with  A(^,y)  for  a 
constant  value  of  X  ^  Xc  •  The  identity  must  then  be  modified  to  read 

y  2  f (y)  =  y2  f(-) 

VI  CONCLUSIONS 

A  method  of  calculating  coefficients  in  L.T.  and  H.T.  series, 
which  makes  use  of  symmetries  and  special  properties  of  the  coefficients, 
has  been  developed  and  applied  to  lattices  of  1,  2,  3  and  4  layers.  The 
method  is  of  most  value  for  homogeneous  lattices  with  small  coordination 
numbers,  but  it  is  useful  for  all  lattices. 

Using  the  series  thus  derived  for  the  one  layer  lattice  and  the 
exact  formula  for  spontaneous  magnetization,  we  have  made  use  of  the 
correspondence  between  the  Ising  ferromagnet  and  a  model  of  sitewise  mono- 
layer  adsorption  with  attractive  lateral  interactions  to  compute  accurate 
isotherms  for  the  latter.  Having  determined  the  consequences  of  this  model 
without  recourse  to  approximate  methods,  which  often  yield  wrong  predictions 
in  the  critical  region,  a  comparison  with  experiment  was  undertaken.  The 
values  of  surface  binding  energy  and  interaction  energy  obtained  from  the 
observed  values  of  critical  temperature  and  pressure  are  the  correct  order 
of  magnitude.  Exact  comparisons  are  not  possible  since  the  experimenters 


. 


1 
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(6,7,8]  did  not  measure  heats  of  adsorption.  However  the  surface  binding 
energy,  Eg,  Is  approximately  constant  for  a  given  surface,  as  one  would 
expect . 

A  comparison  of  theoretical  with  experimental  isotherms  for 
krypton  on  sodium  bromide  shows  qualitative  agreement.  However  the 
experimental  variation  of  transition  pressure  with  temperature  is  less 
than  the  theoretical,  and  the  experimental  magnitude  of  the  discontinuity 
in  6  is  less  than  the  theoretical.  The  former  discrepancy  is  probably 
due  to  neglect  of  vibrational  degrees  of  freedom  of  the  adsorbed  atoms  and 
the  latter  to  imperfections  in  the  substrate  lattice.  The  agreement  is 
reasonable  considering  the  simplifications  of  the  model. 

Calculations  for  the  n-layer  Ising  model  are  complicated  by  the 
fact  that  the  lattice  is  inhomogeneous  for  2  <  n  <  oo  .  The  approximate 
variation  of  critical  temperature  with  number  of  layers  is  qualitatively 
similar  to  that  observed  in  certain  experiments  [11,12,13],  and  a  precise 
correspondence  is  not  to  be  expected  for  this  simple  model.  For  the  two 
layer  lattice  it  was  possible  to  estimate  the  form  of  the  singularity  in 
susceptibility.  It  appears  quite  possible  that  the  two  layer  lattice  has 
the  same  type  of  singularity  as  does  the  one  layer  lattice,  (l-T/Tc)  » 

rather  than  a  form  intermediate  between  one  and  infinitely  many  layers. 

This  result  is  in  agreement  with  Domb's  view  that  the  nature  of 
such  singularities  depends  upon  the  dimensionality  of  the  lattice.  It 
would  be  interesting  to  examine  the  singularities  in  spontaneous  magnet¬ 
ization  and  specific  heat  but  this  is  not  possible  with  our  data. 
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Appendix  I  -  Functions  of  Power  Seriea 

Consider  the  functions 


Fix) 


00 


L _ i 

m=l 


F 

m 


m 

x 


and 


g(*l  =  X  an  x' 
n=l 


We  wish  to  find  a  power  series  expression  for  F(g(x)).  By  substitution 


00  00 

F(glx»  =  Fo  +  £  Fn(  J  ..." 

m=  1  n=  1 


=  F  + 
o 


=  F  + 
o 


oo  oo 

X  F-  (  X  ar(m)  '  ) 

m=l  r-m 


(I.D 


X 


r=l  m-1 


X 


is  the  desired  series.  This  is  a  feasible  method  of  obtaining  the  series 
since  the  sum  in  brackets  is  finite. 


The  numbers  a  'm'  s  defined  by 


00  00 


r=m  r=l 


have  the  following  properties: 

a  («)  _  f 1  .(-) 

r  /_  r-i  i 

1=1 

(m )  _  . _  , 

av  '  =  0  if  r  <  m 

r 


(1.2) 


(1.3) 
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Equation  (1.2)  makes  it  possible  to  express  a^”^  in  terms  the  set  (an) 
for  any  particular  value  of  r  . 

The  particular  case  used  in  this  thesis  is 


in 


00  00 

(1  +  I  an  -  I  br  X"  ■ 


n=l  r=l 

and  by  the  above  equations  we  have 


m=l 


(-l)m+1  a  (m) 

-  r 

m 


This  yields: 


(1.4) 


b_  = 


b_  = 


b,  = 


b.  = 


b„  = 


1  2 

a.  -  —  a, 

2  2  1 

1  3 

a  -  a  a.  +  -  a' 

3  2  i  3  i 

12  2  14 

^  a0ax  • 

a5  '  Vl  -  a3a2  +  a3a?  +  a2al  '  a2al  +  5  “l 

12  2  13 
a,  -  a  a  -  a.  a  -  —  a  +  a,  a  +  2a  a  a  +  —  a^ 
o  4  2  23  4  1  321  32 

3  3  2  2  4  1  6 

~a3al  “  2  a2al  +  a2al  "  6  ai 


(1.5) 


The  equations  derived  this  way  can  be  checked  by  choosing  an  example  like 

CO 

/ ,  \-n  V"1  f  n+k“l  \  k 

(i-x)  =  2,  (  k  ;  x  ■ 
k=o 


W 

in  (l”x)  n  =  n  ^ 
k=l 


x 
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Appendix  II  -  The  Expansion  Coefficients 

We  wish  to  evaluate  the  expansion  coefficients  in 


(&1  ■  l  ’’  - 


By  differentiation  and  a  little  manipulation  we  obtain 

00 

-  2n  =  (1-x2)  £  C/  r  x*'1  .  (II. S) 

r=l 

which  leads  to  the  recursion  formula 

(r+1)  C®+1  =  (r-l)  c”m1  -  2nC"  ,  r  >  0  (II. 3) 


with  initial  values  C™  =  1  ,  =  -2n  . 

Then  by  multiplying  (II. 2)  by  [(l-x)  /  (l+x)]— ^  we  obtain,  using  (II. 3), 
two  more  formulas: 


n  Cn+1 
r 

*  (n+r)  C*  -  (r-l) 

(II. 4) 

„  rn-l 
n  C 

r 

=  (n-r)  C*  -  (r-l) 

(II. 5) 

Subtracting  these 

last  two  yields 

«  rn+1 
n  C 

r 

_  ^n  _n-l 

=  2r  C  +  n  C 
r  r 

(II. 6) 

The  initial  values 

=  0  and  =  2(-l)r  for  r  >  0 

r  r 

are  easily  deduced. 

The  coefficients  were  evaluated  numerically  from  (II.6)  using  an  IBM 
1620  computor. 
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Appendix 


in  A 


in  n 


in  A 


III  -  Series  Expansions 

One  Layer 

4 

=  y  x 

+  y2(2x6-|  x8) 

+  x12) 

+  yU(x8+l8x10-85x12+H8x14-^  X16) 

5,o  10  12  14  16  Qr7.  18  1476  20 v 
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+  u^(  18v2-56v\~^  v6) 

+  uU(  1+50v2-317v\-526v6--‘™;  v8) 

+  u^(138v2-1512v\-4484v6-5072v8+^~  v*°) 

+  u6(2+370v2-6495v^+30,248v6»58>238v8+49,992v10-^a~^  v12) 

5 

+  u7(  986v2-25 , 904v^+ 176,008v6- 509,  520v8+724, 564v10-501 , 232v12+^-^8~v1^) 

+  u8(|+2,586v2-97,789v^+922,478v6-^~i^  v8+7,806, I06v10-8,791, I59v12 

.  14  4,768,955  16, 

+5,093, 150v  — v  ) 


Two  Layers 
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in  a  =  u(|v2)+  U2(l0v23f  v**) 
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Three  Layers 
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Four  Layers 


,  A  /I  5  1  6, 
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Appendix  IV  -  Transition  Pressure  for  Critical  Adsorption. 

The  variation  of  transition  pressure  with  temperature  is  given 
by  the  Clausius-Clapeyron  equation, [4? 3 » 


dp 


o 

dT 


L 

T(V  -V  ) 
g  a" 


(IV. 1) 


where  V  and  V  are  the  volumes  of  gaseous  and  adsorbed  phases,  and 
g  a 

L  =  T(S  -S  )  is  the  latent  heat  of  the  transition.  We  may  neglect  V 
'  g  a  a 


and  treat  the  gaseous  phase  as  ideal. 
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i .  e. 


Then 


or 


a 

v  -  51 
g  Po 


dPrt  lp« 

o  _  o 


dT 

dP, 


RT, 


LdT 


RT 


(IV. 2) 


(IV. 3) 


Since  the  enthalpy  function,  H,  has  the  property 
dH  =  TdS  +  Vdp  , 

then  for  a  transition  at  constant  T  and  P  we  have 


L  =  H  -  H 
g  a 

Using  the  thermodynamic  relation 


we  obtain 


TdS  =  CpdT  -  T  [jfcj  dp 


H ■  I  VT  +  1 1  V'T  (!)  }  «  dT 


+  Ho  • 


(iv. 4) 


Using  (IV.2)  the  second  integral  vanishes  for  both  phases.  Thus 

rT 

L  =  /  (C  -  C  )  dT  +  Lo 

J  PS  pa7 
o 


(IV. 5) 


Assuming  for  simplicity  that  both  phases  may  be  treated  classically,  the 
specific  heats  will  be  constants  and 

L  =  (C  -  C  )T  +  Lo  . 

Pg  Pa7 


Integrating  (IV. 3)  then  yields 
in  p 


(C  -C  ) 

PSD  jal  in  T 

R 


Lo 

RT 


+  const 


Taking  anti logarithms  yields  eq.  (36)  of  the  text. 
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